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ABSTRACT 

Recent observations of entropy of intracluster medium (ICM) at high redshift are 
not consistent with a redshift independent entropy floor. In order to interpret these 
observations, we study models in which the entropy of ICM is enhanced at the epoch 
of cluster formation (2/), which then passively evolves as a result of radiative cooling 
until the epoch of observation {zq)- We confirm that a z-independent entropy floor is 
incompatible with the observations. We find that models in which energy deposition 
into the ICM increases with redshift (scaling approximately as ~ M^(^{1 + z/)^-^) are 
consistent with observations of low and high redshift clusters. Possible sources of such 
non-gravitational heating are briefly discussed. 
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1 INTRODUCTION 

X-ray observations of ICM in galaxy clusters have shown 
its entropy to be larger than that expected from pure 
gravitational processes (Ponman et al. 2003, and references 
therein). These observations have led to the suggestion that 
the entropy of ICM has been enhanced. Possible mechanisms 
for this include 'preheating' the gas before its infall into the 
cluster potential (Kaiser 1991), shock heating during the 
accretion into the cluster (Tozzi & Norman 2001; Babul et 
al. 2002; Voit et al. 2003), radiative cooling (Voit & Bryan 
2001), and energy input from supernova driven winds (Wu 
et al. 2000) or active galaxies (Velageas & Silk 1999; Nath 
& Roychowdhury 2002). Recent observations of lack of an 
entropy core in the ICM in low mass clusters have been how- 
ever difficult to interpret in terms of an entropy floor that 
is expected in 'preheating' scenarios (Ponman et al. 2003; 
Pratt & Arnaud 2003). 

Very recently, observations of the ICM entropy at high 
redshift have shown that a z-independent entropy floor is dif- 
ficult to reconcile with observations (Ettori et al. 2004) . The 
observed entropy (defined as S = T/nJ ) is found to scale 



(where E^ = H{z)/Ho 



with redshift as SE^'^ ex (1 -|- 
is the ratio of Hubble constants at redshifts z and zero). 
In the ACDM universe with Slo = 0.3 and Q,h = 0.7, 
Ez ex (1-1- z;)"'^, which would have led to a stronger scaling 
of entropy (SEj ) with redshift than observed if entropy is 
determined by a z-independent floor. 

It is important to compare these observations with ex- 
pectations from simple considerations of redshift scaling of 
various parameters. For example, in the absence of entropy 
enhancement due to any of the processes mentioned ear- 



lier, one expects the clusters forming at high redshift to be 
denser than their low redshift counterparts, thereby making 
the entropy lower than at low redshift for a given cluster 
mass. It is not clear how the requirement of extra entropy 
should scale with redshift, and this is what we propose to 
discuss in the present paper. 

One important point with regard to the redshift scal- 
ing is that these scalings refer to the redshift of observation 
of a cluster, Zo, whereas any physical model to understand 
the scaUng must involve the formation redshift {z[ > Zo) 
of the cluster, as well as the evolution thereafter until Zo 
(e.g., Kitayama & Suto 1996). In the hierarchical structure 
formation theory, haloes merge to form bigger haloes as 
time elapses. Although haloes also become bigger by slow 
accretion of matter, merging between haloes is the domi- 
nant mechanism of driving haloes upward in the hierarchy 
of structure. In this scenario, there are epochs in between 
major mergers during which haloes evolve slowly by accre- 
tion at a low rate, during which the gravitational potential 
does not change substantially. In other words, when a clus- 
ter is observed at Zo it is reasonable to assume that it has 
been evolving with a reasonably stable gravitational poten- 
tial since an earlier epoch, which can be termed its formation 
redshift zj > Zo- Observational scalings of entropy with Zo 
must therefore be translated to requirements of excess en- 
tropy which is likely to depend on the formation redshift Zf, 
and the duration of evolution between Zf and Zo- 

In this Letter, we build models of galaxy clusters formed 
at high redshift, based on simple extrapolation of those used 
for clusters formed at the present epoch. We then study the 
passive evolution of ICM as a result of radiative cooling be- 
tween the epochs of formation observation. We enhance the 



© 0000 RAS 



Biman B. Nath 



entropy of ICM in various ways at its formation redshift and 
compare with the observations. We then discuss the plau- 
sible models of entropy enhancement which are consistent 
with observations. 

Throughout the paper we use the ACDM cosmological 
model with f^A = 0.7, Q.m = 0.3 and h = 0.7. 



2 ICM AT HIGH REDSHIFT 

We study the distribution of intracluster gas assuming that 
the total gas mass is very small compared to the total mass 
of the cluster, and that the gravitational field is dominated 
by the dark matter. We first determine a 'default profile' 
that the ICM is expected to achieve, in the absence of any 
non-gravitational modifications to its entropy (§2.1). 



2.1 Gas profile without entropy eniiancement 

We assume the ICM to be initially in hydrostatic equilib- 
rium with the gravitational potential of the dark matter. 
The dark matter is assumed to follow the universal Navarro- 
Frenk- White (NFW) profile, pdm = Ps/((r/r,)(l + (r/r,))), 
where rs = r-vir/c, 'c' being the concentration parameter and 
Tvir being the virial radius (Navarro et al. 1997). The virial 
radius is fixed by the overdensity estimated from spherical 
collapse model, for the appropriate redshift of the formation 
of the cluster a/. The overdensity is approximately given 
in a ACDM universe by A{z) = IStt^ -|- 82a:: - 39a;^ , where 
X — ^rn{zf) — 1 (Bryan & Norman 1998), and whose value 
at z = is approximately ~ 100 {flm(z — 0) being assumed 
to be 0.3). The characteristic density ps for a cluster with 
virial mass Md is then given by. 



pB = c 



Mc. 



47fj,3 



ln(l + c) - 



(1 + c) 



(1) 



We use the results of the simulations by Bullock et al. (2001) 
to estimate the concentration parameter as, 



c = 9 



M,. 



1.5 X 101^/i-iMg 



(1 + ^/) 



(2) 



To facilitate the comparison of our results with observa- 
tions which are often expressed in terms of the radius i?2oo 
where the overdensity is 200, we compute this radius and 
present our results in its terms. 

We assume, in accordance with the results from nu- 
merical simulations by Loken et al. (2002) (see also Bor- 
gani et al. (2004); Yepes et al. 2004), that gravitational in- 
teraction of baryons with the dark matter gives rise to a 
'universal temperature profile', which is approximated by 
T{r) = 1.33Tew (1 + 1.5r/r^i^)"^ -^ where Te„ is the emis- 
sion weighted temperature of the cluster. We estimate the 
emission weighted temperature from the empirical mass- 
temperature relation for local clusters and include a red- 
shift dependent factor in this relation that is expected from 
theoretical considerations. The empirical relation found by 
Finoguenov et al. (2001) is given by, 

1.58 
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(3) 



IkeV/ 

where M500 is the mass contained within the radius where 
the overdensity is 500. This relation is valid for clusters with 



M500 > 5 X 10^ M0 and is a reasonable assumption for the 
default profile for all clusters. 

We have used the temperature profile of Loken et 
al. (2002) in our calculations because of the convenient ana- 
lytical fit. It is possible that profiles obtained by other adia- 
batic simulations (e.g., Borgani et al. (2004)) are somewhat 
different than the one used here. A temperature profile which 
is shallower in the inner region than that used here may lead 
to a slight overestimation of the required excess entropy in 
the final result (for details of the implications of the 'univer- 
sal temperature profile', see Roychowdhury & Nath 2003). 

The expected scaling of the mass-temperature relation 
with redshift Zf of formation of clusters is T oc E^i M^'"^, 
where E,^ = H{zf)/H{z = 0) = [no(l + Zff + Aol^'''' (e.g, 
Afshordi & Cen 2002). We therefore multiply the observed 
mass-temperature relation (equation 3) with a factor of E^^ 
on the right hand side to determine the default profile at the 
formation redshift. 

The initial gas density profile is then assumed to be 
determined by the hydrostatic equilibrium condition, -£: = 
—pgGM{r)/r^, where the temperature profile is determined 
as explained earlier. We normalize the gas density distribu- 
tion in such a way that the total gas mass within the radius 
-R200 is a constant fraction /g,200 = 0.133 of the total cluster 
mass, as recently found by Ettori (2003) for high and low 
redshift clusters. 

We express our results in terms of the emission- weighted 
temperature of the ICM, calculated with the help of the 
Raymond-Smith code for a metallicity oi Z = O.SZq. 



2.2 Enhancement of entropy 

We assume two simple methods of entropy enhancement to 
the initial gas profile. Let us define a = P/p'' as the 'entropy 
index', where 7 is the adiabatic index. The entropy per unit 
particle can be written in terms of this entropy index as 
s — constant + (7 — l)fcs ln((j), where fc_g is the Boltzmann 
constant. In terms of this entropy index, we assume that in 
the first type of entropy enhancement. 



a{r) = iaax{ao,(7i{r)) , 



(4) 



where (7i{r) refers to the initial gas profile and ao is a con- 
stant, for all radii and for all clusters. Physically, this corre- 
sponds to an entropy fioor which may be caused by preheat- 
ing. Below we will refer to this model as the 'preheating' 
case. 

In the second type of entropy enhancement, we assume 
that, 



a{r) ^riai{r), 



(5) 



where r/ is assumed to be a constant for the given cluster, 
for simplicity. Physically this type of entropy enhancement 
would correspond to energy deposition into the ICM after 
its infall within the cluster potential well. The excess energy 
per particle thus deposited, Ae = TAs oc T— oc T(?7 — 
1). Since gas temperature is a weak function of radius, this 
roughly corresponds to a near uniform deposition of energy 
throughout the cluster, for a constant value of 77. We will 
refer to this as the 'post-infall heating' case in the following. 
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Zo)"^'^(Mci/lO"M0)~''°^ which is accurate within 5% for 

Zo< 1. 



1.2 1.4 1.6 1.6 2 

(I+z„) 

Figure 1. The duration of evolution, between the redshift of for- 
mation, Zf and that of observation Zo of clusters of mass 10^'* Mq 
is shown against 1 + Zo- The thick solid line corresponds to the 



median value of distribution of , 



and the shaded region corre- 



sponds to 30% and 70% values of Zf . The dashed line corresponds 
to the median value of Zf for clusters of mass 10^"" M©. 



2.3 Halo formation epoch 

Since haloes increase their mass by mergers and steady ac- 
cretion, it is difficult to assign a single epoch of formation 
for a given mass. Lacey & Cole (1993) defined the halo for- 
mation epoch Zf for a halo of mass M to be the epoch at 
which the halo has a mass greater than M/2 for the first 
time. After this epoch ('of formation') the halo is supposed 
to attain the total mass M by steady accretion of small ob- 
jects. Lacey & Cole (1993) then determined the distribution 
function for the probability that a halo of mass M that is 
observed at z was formed at Zf, or had a mass greater than 
M/2 for the first time at Zf. 

For our calculations, we define the halo formation epoch 
as the epoch when the halo has a mass greater than |M for 
the first time. This is motivated by the results of numerical 
simulations which show that gravitational potential does not 
change much after 75% of the total mass is in place (Navarro 
et al. 1995; see also Babul et al. 2002). 

For our calculations, we use the power spectrum for 
ACDM cosmology normalized by ag = 0.9, and for growth 
of linear perturbation in this cosmology, we use the fit by 
Lahav et al. (1991).. Using the probability distribution we 
can determine different percentiles of the formation epoch 
Zf. We show in Figure 1 the expected time of evolution, 
i.e. the time spent between the formation redshift Zf and 
when it is observed at Zo, as & function of Zo for clusters 
of 10^* Mq. The solid line represents the median value of 
the distribution in zf and the shaded region the 30% and 
70% percentile values. The dashed line shows the median 
value for Zf distribution for clusters of mass 10^^ M0. It is 
clear that evolution is more important for clusters observed 
at lower redshifts. 

We find that a convenient fit for this duration is 
given by (for the median value) At ~ 2.5 x lO'^yr (1 + 



3 EVOLUTION OF GAS IN CLUSTERS 

We divide the gaseous sphere (which resides within the po- 
tential well provided by the dark matter, Mdm ~ Mt, which 
is the total matter) into mass shells which are designated by 
gas mass contained within them, Mg{r). In terms of the en- 
tropy index introduced earlier, the hydrostatic equilibrium 
condition reads. 



dP 

dMg 

dr 

Ufa 






ATvr 



liPM 



-1/7 



(6) 



The default profile is computed using these equations and 
the universal temperature profile. We set the boundary con- 
dition that the gas pressure at the mass shell which was 
at a radius R200 for the default profile, and for which 
Mg/Mci = /a, 200 = 0.133, remains a constant. 

The entropy of the gas is then enhanced instantaneously 
(according to equations 4 and 5). This enhancement changes 
the density and temperature profile of the gas (pushing it 
outward to some extent; see below). 

We assume that after the initial (impulsive) enhance- 
ment, intracluster gas cools radiatively and evolves in a 
quasi-hydrostatic manner. For subsequent evolution as a re- 
sult of radiative cooling, if A(f ) denotes the cooling function, 
the specific entropy evolves with time as (see also Kaiser & 
Binney 2003), 



ds 
dt 



7 



1 , 1 dcr 
— 1 G dt 



A(r), 
r 



(7) 



where the electron density Ue = 0.85p(r)/7np. We use the 
cooling function appropriate for metallicity Z — 0.3 Zq from 
Sutherland & Dopita (1993) for which Nath (2002) provided 
a fit. The equation 7 is integrated for each mass shell for a 
small time interval St and the entropy index added to this 
shell after the time interval is. 



Act = -(7- 1 



AT , 
keT 



(8) 



The equations of hydrostatic equilibrium are used again to 
determine the density and temperature profile at time t + St. 



4 RESULTS 

We first show the results for the case of preheating, i.e. 
the case of an entropy fioor (equation 4), in Figure 2. It 
shows the observed entropy at 0.1i?2oo for clusters with 
different emission weighted temperatures and observed at 
Zo = 0.0,0.5, 1.0 for an assumed value of ctq = 0.27 x lO^'* 
erg cm^ g~^ • (For comparison, Tozzi & Norman (2001) 
used a value of 0.3 x lO^** erg cm^ g~^ for their calcula- 



2/3 



247.9 keV cm^ 



tions.) This corresponds to S* = T/i 
(assuming rie ~ 0.875p/mp). It is constructed out of data 
points calculated for clusters of different masses observed at 
different epochs. The extreme left points in each set corre- 
spond to clusters of 5 x 10^'^ M© and the extreme right points 
use clusters with 5 x lO^** Mr. The thick solid line uses the 
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z^=(r.3(In)" 



0,2 0.4 0.6 0.8 1 1.2 

Log T,, 

Figure 2. Plot of entropy (at 0.1i?20o) for clusters with different 
Teiij and observed at redshifts Zo = 0., 0.5, 1.0. Tlie tliick solid 
lines in each redshift set correspond to the median value of dis- 
tribution of Zf and the shaded region (with dots) to the range 
between 30% and 70% values of Zf. For comparison, the initial 
entropy of the same clusters are shown with dashed lines (and la- 
beled 'in'). The three slanted dotted lines show the observed scal- 
ing of entropy-temperature relation, Ej S oc TO. 65^]^ _^ ^^Y^.:! ^ 
corresponding to the three different redshifts Zo = 0., 0.5, 1.0. 



median value of distribution of zj, and the shaded region 
corresponds to the range between 30% and 70% values of 

Figure 2 leads to a few important conclusions. Firstly, 
it is clear that a z-independent entropy floor is inconsistent 
with both the observed S-T relation at low redshift and its 
scaling with redshift. The S-T relation that follows from it 
is in general shallower than the observed S oc T^-^^ scaling 
at the present epoch (Ponman et al. 2003) . This result has 
been earlier noted by Borgani et al. (2001) from numerical 
simulations. The evolution of ICM between the formation 
redshift and that of observation makes this problem worse. 
This is because gas in massive clusters (with higher Tew) 
loses entropy at a faster rate. This is elaborated in Figure 
3, which shows a few examples of evolutionary tracks of en- 
tropy with redshift for different clusters with various Zf . The 
solid and the dotted lines show the evolution of a massive 
and poor cluster, respectively, with the same Zf . Clearly the 
ICM in massive cluster loses entropy faster than the ICM 
in low mass clusters, making the observed S-T relation in 
Figure 2 flatter than the initial S-T relation at formation. 
Moreover, since the duration of evolution (between Zf and 
Zo is longer at low redshifts, the flattening of S-T relation is 
prominent at lower redshifts. 

Secondly, the evolution of the S-T relation with red- 
shift is stronger than observed (shown with dotted lines). 
This has already been pointed out by Ettori et al. (2004). 
Our calculations show that the problem is made worse than 
previously thought as a result of evolution of ICM between 
Zf and Zo- 

We then turn to the case of post-infall heating. We have 
found that assuming 7] to be independent of cluster mass and 
redshift produces a steeper S-T relation than observed, and 



7 xlO" Ma.Z| = 0.5 

10" M3,.z,-0.5 

10" M,.z,-1.0 

J I , , , I ^ 

1.4 1.6 



(l+z) 

Figure 3. Evolution of entropy at O.IR200 for clusters of different 
masses with different redshifts of formation is shown. Dotted and 
dashed lines correspond to clusters of mass 10^'* Mq with forma- 
tion redshift 0.5 and 1.0 respectively, and solid line correspond 
to, 5 X lO^"* Mq with Zf = 0.5. 




Log T, 

Figure 4. Same as in Figure 2 but for the case of 'post-infall' 
entropy enhancement, with rj oc M~, ' (1 -(- Zf)'^''^ . 



a scaling with redshift slower than observed. We therefore 
used simple power-law scalings for rj with cluster mass and 
redshift in the form, ri{Mci,Zf) oc M"i{l + Zf)" , and de- 
termined the approximate values of a and /3 that produce 
results consistent with observations. 

We show in Figure 4 the results for the 'post-infall' heat- 
ing case, with »? = 10(Mc;/5x10"Mq)"°-^(1-|-2/)2-2, along 
with the observed S-T relation and its scaling with redshift 
(shown with dotted lines). We flnd that this form of de- 
pendence of 77 with cluster mass and formation redshift is 
reasonably consistent with observations. 

Finally, we show a few examples of the entropy pro- 
files, before and after the entropy enhancement and after 
the passive evolution until the epoch of observations, for 
a couple of clusters with varying masses in Figure 5. The 
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0, z,=0.16 

default (z = Zf)" 
heated (z = Zf) 
evolved (z = z^ 

1 -0.5 

Log r, 

2/3 

Figure 5. Entropy (S = T/nJ ) profiles are shown at differ- 
ent epochs for two types of entropy injection and for two cluster 
masses. The left panels represent clusters of mass 5 X 10^'' Mq 
with Zf = 0.2 and the right panels, 7 X 10^'* Mq zj = 0.16. The 
upper panels show the case of z-independent entropy fioor (equa- 
tion 4) and the lower panels show the case of 'post-infall' heating 
(equation 5). The solid curves show the initial profile, the dashed 
lines show the profile immediately after entropy injection and the 
long-dashed curves show the result of evolution until the redshift 
of observation at Zo =0.. 



solid line shows the initial entropy {S = T jnj ) profile at 
formation epoch. The dashed curve then shows the profile 
after entropy enhancement. The upper panels show the re- 
sult of using equation ( 4) for entropy injection and the lower 
panels use equation (5). As expected, the intracluster gas is 
pushed out as a result of this entropy injection. Finally the 
long dashed curves show the result of evolution until the 
redshift of observation (zo = in this case) . The left panels 
show the case of M = 5 x 10^^ Mq and the right panels that 
of M = 7 X 10" Mq. The redshift of formation for these 
two cases have been chosen to be the corresponding median 
values of distribution of zf (0.16 and 0.20). 



5 DISCUSSION 

We discuss the implication of the result arrived at the last 
section that energy deposition in ICM with r\ oc M^°'^(l -|- 
2:/)^'^ is consistent with observations of low and high- 2: clus- 
ters. Firstly, since the excess energy deposited per particle 
is proportional to T(r] — 1) ~ Tr\ (since 77 > 1), our assump- 
tion of a constant r\ for a given cluster roughly amounts to 
uniform deposition of energy throughout the cluster, since 
temperature is a weak function of radius (it decreases by 
much less than an order of magnitude from the core to the 

outermost radius for the default profile) . Secondly, since T oc 

2/3 
M^/ , the variation of r\ with Md implies that the required 

1 /3 

energy per particle Tr\ ~ 10 keVM^/ . Since gas mass is pro- 
portional to the total mass, the above result implies that the 
total required amount of energy Ai? oc AeMci oc M\f^ . 

The required scaling of r] with redshift (oc(l-|-2;/)^'^) 
is mainly due to the fact that densities are larger at higher 
redshift, reducing the default entropy level, and thereby in- 



creasing the requirement. This increase is somewhat com- 
pensated by the relatively short duration of evolution (be- 
tween Zj and io) at high redshift ( see Figure 1) but not 
significantly. 

It is of course possible that the enhancement of entropy 
is more complicated than represented by equations 4 and 5. 
Especially the form of energy deposition into the ICM after 
its infall can either depend on the radius and/or time, in 
which case the required scaling of r\ would be different from 
the results presented here. A non-trivial radial dependence 
of r\ may change the entropy profile substantially and change 
the form of passive evolution thereafter. For example, in the 
case of a strong radial dependence, the entropy in the inner 
region may increase substantially, with the possibility of a 
negative entropy gradient ushering convective fiows (e.g., 
Roychowdhury et al. 2004). 

Our assumption of an instantaneous entropy enhance- 
ment followed by passive evolution has been admittedly sim- 
ple, since our goal has been to interpret the general depen- 
dence of observed entropy with redshift. It is however pos- 
sible that the period of entropy enhancement is substantial 
(see, e.g., Roychowdhury et al. 2004), but modelling such 
cases introduces another free parameter (of enhancement 
timescale) in the calculation. Although detailed modeling of 
such complicated forms of heating may change the scaling 
power index of r] with redshift to some extent, but is not 
expected to substantially change the strong redshift depen- 
dence of r\, which mainly arises from densities being larger 
at high redshift. 

It is interesting to note that the final entropy profiles 
in the case of near-uniform heating after infall (see lower 
panels of Figure 5) are consistent with recent observations 
of positive gradient of entropy even at small radii for poor 
clusters, instead of a core as expected from the preheating 
scenarios (upper panels of Figure 5) (Ponman et al. 2003; 
Pratt & Arnaud 2003). The near isentropic core in the case 
of an entropy floor is consistent with the numerical simula- 
tions of Borgani et al. (2001) (considering that the value of 
the entropy floor used here is larger than assumed by them), 
and also Babul et al. 2002 (who used a value of ao that is 1.7 
times that used here). We must however remember that a 
strong radial dependence of heating of ICM may change the 
profile by depositing more energy (and entropy) in the cen- 
tral region as compared to the outer region. In other words, 
the recent observations of positive gradient of entropy at 
inner radii is consistent with a (post-infall) near-uniform 
energy deposition into the ICM throughout the cluster. 

As far as the plausible sources of heating which can ex- 
plain the observations of high redshift ICM, it is interesting 
to note that rate of heating by active galaxies is expected 
to rise with redshift (Nath & Roychowdhury 2002) mainly 
because of the fact that quasar activity has declined in the 
universe in the recent past. The amount of energy deposited 
is also reasonably consistent with estimates from heating 
from buoyant bubbles from active galaxies (Roychowdhury 
et al. 2004) . In this regard, one should also consider the pos- 
sibility of supernova driven winds as the source of heating, 
since the total energy deposited in this case is also expected 
to increase with the cluster mass (total mass of galaxies, 
from which the winds would emanate, being proportional 
to the total cluster mass). Moreover since the average star 
formation rate in the universe increases with redshift ( for 
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2: < 1), it is possible to achieve a strong redshift dependence 
in this case. We however note that detailed calculations have 
shown supernova driven winds to be inefficient sources of 
heating (Wu et al. 2000). 



Yepos, G., Ascasibar, Y., Sevilla, R., Gottlober, S., & Miillcr, V. 
2004, in Proceedings of lAU Coll. no. 195 { A. Diafcrio, Ed.), 
in press (astro-ph/0404556) 



6 SUMMARY 

We have calculated the expected evolution of entropy of 
ICM (at 0.1i?2oo) with redshift (of observation Zo), taking 
into account passive evolution of ICM between the epochs 
of formation and observation of the cluster. We found that 
the adoption of a z-independent entropy floor exacerbates 
the problem of accounting for the observed evolution of the 
entropy-temperature relation with redshift. If the entropy is 
enhanced instead with the help of a heating source after the 
infall of the ICM, we find that the total energy required to 
be deposited scales as AS oc Mlf^(l + z/)^'^. 

It is a pleasure to thank S. Roychowdhury and S. Srid- 
har for stimulating discussions, and the anonymous referee 
for valuable suggestions. 
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